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baryon number violation in r or higher generation quark decays, at any future machine. 
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I. INTRODUCTION 

Proton decay was suggested 30 years ago as a means to 
probe physics at very high energy scale, such as p| Grand 
Unified Theories (GUT). Baryon number violation is also 
a key ingredient [J to understand the baryon asymmetry 
of our Universe. Experimental search for proton decay 
has yielded 3] only limits so far, albeit very stringent 
ones, and improvement requires huge underground facil- 
ities. 

More recently, in an attempt to solve the hierarchy 
problem £|, a new scenario has been proposed, which 
assumes the existence of large compactified extra dimen- 
sions. It has been observed that, having no test of gravity 
below the millimeter scale, the radius of compactification 
can be much larger than Mp lanck , provided the extra di- 
mensions are accessible only to gravitational interactions. 
Standard Model fields are instead localized on a three- 
dimensional wall whose inverse thickness L _1 can be of 
order TeV. Denoting by V n the volume of the extra space, 
with n number of compact extra dimensions, and by M 
the fundamental scale of gravity, the observed Planck 
mass is given by Afp lanck = M n+2 V n . By asking that 
n > 2, the fundamental scale M could be close to the 
electroweak scale. 

In this new scenario, proton stability could be problem- 
atic. It is not possible anymore to request that the gauge 
bosons which mediate baryon number violating (BNV) 
processes to be super massive |j| . A new mechanism has 
been suggested in Ref. where leptons and quarks are 
localized at two different positions in the extra space. 
The four-dimensional effective coupling of the higher di- 
mension operators responsible for proton decay are then 
suppressed by the small overlap of the quark and lepton 
wave- functions along the extra dimensions. 

It has also been proposed @ that by localizing differ- 
ent generations, and at the same time placing left and 
right components at different positions along the extra 
space, one can obtain in a natural way quark mass ma- 
trices which reproduce a realistic mass spectrum and the 
quark mixing matrix. In Ref. Q it has been shown 
that a split fermion scenario, disfavored by constraints 
on flavor changing neutral currents (FCNC) mediated by 
Kaluza-Klein gauge bosons, is still possible if the three- 



dimensional wall has inverse thickness L _1 > 100 TeV. 

Having these scenarios in mind, it is reasonable to 
ask if BNV processes involving higher generations, or 
right handed components, could be faster than the pro- 
cesses involving the first generation, which we know to 
be strongly suppressed. One can imagine an extra space 
geography where the second and third generation quarks 
have big overlap with the leptons while the first genera- 
tion quarks are still sufficiently far away from the lepton 
localization position. 

It is therefore important to give an estimate of how 
fast BNV processes involving higher generations could 
be, and how left or right handed dynamics may affect 
them. In Ref. it was suggested that one can estimate 
BNV t decay branching ratio by using proton decay con- 
straints. Applying the same approach, we estimate BNV 
processes involving higher generations. 

Let us write down all possible dimension six operators 
that violate baryon number, but respect the symmetries 
of the Standard Model (SM) HQS: 

°abcd = (d% aR Uf3 bR )(% 7C JjdL)e a 7 eij, (1) 

°abcd = (^ a aLl3f3bL)(u C 7cR £ dR ) e^eij, (2) 

(3) C C 

°abcd = {^iaaLljPbL)(.q C k 1 cL^ldL)^af3^i 3 ^kU (3) 
°abcd = (QiaaLQj0bL)(qky c L^ldL)£ a 0-y 

x(re)y-(fe)fe;eifc, (4) 
°abcd = (dl aR U0 bR )(u^ cR £dR)e a 0j, (5) 

°abcd = (^ C a aRU0bB)(d C lc R^dR)e a 7 , (6) 

where a, (3, 7 are color indices, i,j are SU(2) indices, 
a,b,c,d are family indices, and e a p-y, tij are totally an- 
tisymmetric [TT|. From Eqs. 0-@ one can build the 
effective interaction 
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E r {n) n (n) / ? n 

° abed'-' abed' \ ' J 

n=l 

by introducing the coefficients C^ cd of dimension d = 

(n) 

—2. The coefficients C abcd are model dependent quanti- 
ties which are sensitive to high energy scales. We shall 
give a model independent estimate of the Cj^L coeffi- 
cients by using experimental limits on nucleon decays. 




FIG. 1: p — > v t tt + generated by C^ udT . 




FIG. 2: p — > A' e ^ e ^r generated by C uuaT . 

The structure of this paper is as follows. In section II 
we estimate BNV r decay branching ratios by using pro- 
ton decay constraints. In section III we apply the same 
approach to higher generation quarks and consider BNV 
processes mediated by tree diagrams. In section IV we 
show that loop diagrams give much stronger constraints. 
In section V we estimate BNV D and B meson decays. 
Our conclusions follow. 



II. BNV OPERATORS AND r DECAYS 

We start by analyzing BNV processes involving t, 
which are intuitively very intriguing. The fact that the 
proton cannot decay into r but only into v T implies 
that BNV decays involving left-handed r are highly con- 
strained, but for right-handed r, BNV processes could 
be in principle faster. However, following Marciano @, 
by considering processes with virtual r, as for example 
p — > "virtualr" — > v T ir + (Fig. 1), one can also strongly 
constrain r decays involving right-handed r. So far the 
search for r~ — > p7r°, pr?, p7r°7r°, pn !], and 
A(A)7T _ ^4| decay modes has given only upper limits. 

Let us calculate explicitly the process of Fig. 1 when 
the BNV vertex involves right-handed r. We use the chi- 
ral Lagrangian obtained from the operators of Eqs. J3J), 
(J5J and 10. Each operator transforms under SUi(3)x 
SUj?(3) like (3,3), hence the AB — 1 Lagrangian is 

% (c« udT Tr 0£tB^t + C R usr Tr 6^BU f ) tr, (8) 

where £ = exp^M//^) with M and B the pseudo- 
Goldstone boson and octet baryon fields, respectively, 
and C^ udr , C^ust are effective couplings for opera- 
tors in Eqs. ©, © and ©. The matrices O and 
6 project out AS = and 1 terms [H [Oj. The 
strong interaction parameter a p is defined by the rela- 



tion (0|e a/ 9 7 d a .RM / 9.Rti 7 i,|p) = a p Up L , with u p the pro- 
ton spinor. A recent lattice calculation gives a p = 
—0.015 GeV 3 , applicable for proton decay. The BNV ver- 
tex in Fig. 1 corresponds to 0i p C^ udr f^pR from Eq. JSJl, 
with ot p now at r scale. A straightforward calculation for 
the mode p — > v T ir + shown in Fig. 1 gives the decay rate 

with the Fermi coupling constant. From the experi- 
mental limit r(p — > p7r + ) > 25 x 10 30 years, one can now 
constrain C^ udT and obtain 

G* dr < 6.1 X 10" 18 TeV" 2 , (10) 

where we have applied the same value for a p as in proton 
decay. 

Using Eq. JHJ, the general SUl(3)x SUa(3) invariant 
chiral Lagrangian involving the M and B fields [l5j , and 
the r width |3j, we obtain the branching ratios, 

B(r^pn a ) < 1.6xHT 4 |G« dr | 2 , (11) 
B(t - pK°) < 7.5 x 10- 5 \Ctsr\\ (12) 
B(r^ An-) < 1.3xlO- 4 |G« sr | 2 , (13) 

where C^ udr , C^ UST are in TeV -2 units. Substituting 
Eq. (prnf> in Eq. (jflj one obtains 

B(t -» pvr°) < 5.9 x 10~ 39 . (14) 

The other two modes involving a strange quark depend 
on the effective coupling C^ UST , which corresponds to 
the —ia p / f v {C^ UST f%pRK ^) term from Eq. ©, with a p 
again at r scale. Considering the process p — ► K°e + v e v T 
of Fig. 2, and using t > 0.6 • 10 30 years for the semi- 
inclusive process p — > e + + anything, one obtains 

G* ST < 1.8 x 10- 13 TeV" 2 , (15) 

where again we have applied the same value for a p as in 
proton decay. Substituting now Eq. (|15|) in Eq. I|12l) and 
Eq. i|13|) we also obtain 

B(r^pK°) < 1.4 xl0~ 30 , (16) 
B{t -» Xtt") < 2.3 x 10~ 30 . (17) 

One can generate r — > py, £+7 decays from Eq. JSJ by 
radiating a photon off the external lines in the effective 
r — * p, £+ transitions, which are analogous to p — > £ + 7. 
Direct radiation from the transition vertex would involve 
higher dimension operators, hence suppressed. The ra- 
diative process should be down by a power of a cm /ir com- 
pared to the modes we have just discussed. Indeed, we 
find the branching fraction, 

B(t - py) < 2.8 x 10- 7 \C* udT \ 2 £ 1-0 x 10- 41 . (18) 

We see that BNV r decays are extremely suppressed, 
and it is unlikely that we will ever observe them at any 
future machine. We now turn to BNV processes involving 
higher generation quarks. 
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FIG. 3: 



n + TT 7r e + (fj, + ) generated by Ct, 



III. C tt , c)be COUPLING FROM n -> 



In the previous section we showed how Marciano's sug- 
gestion of considering virtual r in proton decays can be 
used to constrain C^ udT and C^ UST . By analogy one can 
then constrain any coupling C a bcd through nucleon de- 
cays, where now the family indices a, 6, c and d can 
belong to any generation. 

As an example, we give in Fig. 3 the tree diagram 
for n — > £ + (pi)ir + (p 2 )tt~ {p3)tt~ (Pi) induced by a C t m 
coupling, where £ + could be e + or fi + . As internal states 
all possible combinations are allowed, where a b (t) quark 
propagator can be replaced by a d or s (u or c) quark. 
The amplitude will then be a combination of terms each 
proportional to a different C a bcd- 

As we are interested in giving only an order of magni- 
tude estimate, we can use t(tl — > e + (/x + ) + anything) > 
0.6(12) x 10 30 years to constrain the C a b c dS, by assuming 
each time that only one single term dominates. Taking 
for example the term proportional to Cum as the domi- 
nant contribution and considering for simplicity only the 

(5) 

effect of O n L, the amplitude can be written as 



M {5) 



VIVA C 



(5) 
ttbl 



2^2 ^ U 
(l + TT-TT+TT-\ [u C in(\ + l5 )b C ] [b R t R i R £ R ] 

+ 7 5 )i c ] ©4(1 - J 5 )d] \n). (19) 



Using the matrix element (0\u c d c d\n) = a r , 



v n , we obtain 



iM^ tmnt = 2V2G F flv: b V t 2 d a n c[ b) 



lib l 



1 



tr(^ 4 ) ^3(1 + 75)^, (20) 



where we have neglected the internal quark momenta pb, 
Pt and p' t (all of 0(1 GeV)) with respect to mb(l GeV) 
and m t (l GeV). 

It is clear from Eq. I|20|l that the process n — > 
£ + Tr + 7r _ 7r _ is highly suppressed, depending on the 6th 
power of Gf and the 18th power of the Cabibbo pa- 
rameter A appearing in the Wolfenstein parametrization 
of the Cabibbo-Kobayashi-Maskawa matrix (CKM). Fur- 
ther suppression comes from four body phase space and 
the hadronic parameter a n . Thus, the experimental con- 
straint r(n — > e + (fi + ) + anything) > 0.6(12) x 10 30 years 
can be satisfied without requiring the effective coupling 



FIG. 4: 



y T Tv° generated by C t tb T - 



Ctm t° be very small. In fact cj. M 
tain for the n — > £ + tt + it~tt~ 



_ can be huge. We ob- 
rate the following estimate 



r£l w [TeV] 



<4xlO- 87 |C t ^| 2 , 



which gives 



4 5 i<10 n TeV- 



(21) 



(22) 



Obviously, the upper limit in Eq. I|22|l is much higher 
than what could be realistic, but it would still seemingly 
suggest that BNV processes like t — > bb~W~£ + could take 
place rather quickly. Assuming now that 0^ M can give 
the dominant contribution, from Eq. I|20ll it is clear that 
the constraint on C tc bt is also very loose, and we obtain 
the rough estimate 



C tM ~C tM ^—<10 7 TeV- 
V cd m t 



(23) 



Decays like t — > cb£ + could then easily compete with SM 
top decays making LHC a BNV machine. 

Before getting carried away, we remind ourselves that 
the above discussion of Fig. 3 is just an example. What 
we need to do is to identify the fastest possible process 
generated by BNV involving higher generation quarks. It 
is clear that we need to reduce the powers of Gf or A. 
We now turn to what we feel is the best process in the 
next section. 



IV. 
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tt(c)be 



COUPLING FROM 



The process n — * £ + it + 7t~7t~ in Fig. 3 is highly sup- 
pressed, because for each tree level virtual W emission, 
there is one power of Gf suppression. It is clear that a 
diagram involving W emission at loop level will be less 
suppressed. In Fig. 4 we give the two-loop radiative cor- 
rections for the process n — > v T ir° involving Cubr] similar 
diagrams describe two-loop radiative corrections for the 
process n — > D e ,^° with r replaced by e and [i. The am- 
plitude for the process n — > D t tt° at tree level is propor- 
tional to C u ddv T ■ When radiative corrections are consid- 
ered, G u ddu T mixes with all the other couplings including 
Cubr- Knowing that C u ddv T is already constrained from 
the tree level calculation, we obtain a rough estimate of 
Cttbr by just isolating the contribution to the amplitude 
given by the diagrams of Fig. 4. We find in this way the 
most stringent constraint on C^ c ^ be . A very similar ap- 
proach has been used in Ref. [lj| to constrain the gauge 



couplings of an SU(2)i x SU (2)2 extension of the Stan- 
dard Model. This theory can in fact give rise to baryon 
and lepton number violating processes via the instanton 
transition associated with the extended gauge group. 

Again we are interested only in order of magnitude 
results, and it is therefore unnecessary to sum up all am- 
plitudes coming from different operators. Let us assume 
(3) 

for example that C\ tbl gives the dominant contribution 

(3) 

by considering only the effect of the operator OLL in- 
volving left-handed components. The amplitude for the 
first diagram in Fig. 4 gives 



M 



(3) 
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.(3) 



S—v* v 2 r 

64 v ub v td ^ttbl 



d 4 h d 4 k 2 



(27r) 4 (27r) 4 



{kl-M^){kl-Ml) 



[J C 7 Q (1 + 7 5 )* C ] [S C 7 /3 (1+ 7 5 )6 C 
[g 7 "(l-7>] [Ht L i c L £ L ] \n). 



(24) 



Using the matrix element (ir°\u c d c d\n) = a n7r w„ |19| . and 
after some simplifications, we get 



2 G%Mw V* b V 2 d a n7r C^ b \ 



d 4 hd 4 k 2 tr [7^7^(1 - 7 5 )] 



(2^(2^)4 (k\-M 2 w )(k 2 -M 2 w ) 

Vnla&p'll^ -75)^ 

(Pb ~ m l)(Pt ~ m t)(p't 2 ~ m t)(pj ~ m T) ' 



(25) 



where k\ and ki are the momenta of the internal W 
bosons, and p b , pt, p\ an d Pi are the momenta of internal 
quarks and lepton. 

Note that our BNV interactions are of four fermion in- 
teraction form, hence not renormalizable. The two loops 
integral in Eq. (|25(l is ultraviolet divergent. To obtain a 
finite result we can introduce a cut-off scale A where we 
assume our effective theory will break down. Setting all 
the external momenta to zero we obtain 



(3) — 2 

where C\ tb \ is measured in TeV . Using the exper- 
imental lower limit on the n — > vtt° partial lifetime, 
T n ^ P7r o > 112 x 10 30 years, we obtain 



< lCT 14 TeV^ 2 . 



(29) 



One can follow the same approach for the operator 
Otm w hich involves only right-handed fields [2(|, allow- 
ing for a weaker constraint on CjSt 

C^l < 1CT 11 TeV" 2 . 



(30) 



Assuming now that 0^ bl gives the dominant contribution 
and using Eqs. (|29|l and H3U|) we obtain 

C& ~ «^^4 S ID" 16 TeV" 2 , (31) 



Vcd I(m b ,m c 



with I(m a , m b ) defined by I a ^(ma 7 m b ) = g al3 I(m a , m b ), 
and 

^-^^^T^^^lO^TeV- 2 , (32) 



where 



Is 



d 4 k 



(2tt) 4 (k 2 - m 2 a )(k 2 - m 2 ){k 2 - M 2 ,) 



(33) 



From Eq. (|32|) we expect t — > cb£ + can have at most a 
branching ratio 



B(t^cb£+)~lO- 5 \C^f<lQ 



-27 



(34) 



(51 

where C tih \ is in TeV . This excludes any possibility 
for discovering BNV top decays at future colliders. 



M n ->Pw° = -2 G 2 F Mw V* b V 2 d aw C^ b \ 

/( l)^ /( 2)™ tr [ 7 a 7r7e7M(1 „ 7 5 } ] 

where /M"^ is defined as 



j(i)af3 



d 4 k, 



h a h p 



(2tt) 4 (k 2 - m 2 a )(k 2 - m 2 )(k 2 



M 2 W 



(26) 



(27) 



with (m a ,rn b ) respectively (pit, mi) and (m b ,mt). The 
integral /W"' 3 has logarithmic dependence in the cut-off 
scale A plus corrections of order A -2 0. Setting the 

scale A = 10 TeV, and considering only the contribu- 

(3) 

tion to the amplitude coming from O ttbi , we obtain the 
following rough estimate for the n — > i>ir° rate 



,(3) 



JTeV] <10- 38 |C, 



381/^(3) |2 



ttbi\ 



(28) 



V. D AND B DECAYS 

In this section we will give rough estimates for BNV 
D and B meson decay branching ratios. Following the 
arguments of the previous sections for BNV r and top 
quark decays we expect BNV D and B processes to be 
also very suppressed. Nevertheless we will still try to give 
some directions for a possible search. 

We start by considering D decays. We first note that, 
to allow for a baryon in the final state, r lepton is forbid- 
den, hence we consider only c — > £ + su, i + du processes 

with £ = e,ii. These lead to the decays D+ -> M+ (A*V 
has more background) and D° — > Y,+£ + , p£ + . The D + 
mode has some advantage because it has a longer lifetime 
and cleaner signature. Since we are not concerned with 
hadronic effects, we emulate the Lagrangian for r decay, 
Eq. JSJ , which is itself inspired by the chiral Lagrangian 
for proton decay. Assuming that hadronic factors are of 
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the same order of magnitude, we find the Lagrangian for 
D meson decays of interest, 

i^C* M D°(f R l R ) + 

JD 

ij^Cucsi (d q (£+ c £x) + ^-D+{A R £ R ) j . (35) 



From this we obtain the branching ratios, 



B(D+ -> A£+) < 3.9 x 10- 4 |C, 
B(D°^W£+) < 3.1 x lO- 4 \C« se \ 
B(D° ->p£+) < 1.1 x IQ- 4 \C, 



R |2 
2 



ucdA 



(36) 
(37) 
(38) 



where the C's are in units of TeV 2 , and we have set 
p . These single baryon plus lepton final states 



ao < a 



are rather distinct, but have so far never been searched 
for. In fact, baryon number conservation instructs one 
to consider baryon pair in final state, which is kinemat- 
ically forbidden. Thus, baryons in D decay final states 
have commonly been ignored. D decay to final states 
containing a baryon, however, only probes baryon num- 
ber violation in second and first generations. 

We do not consider D s decays, as it offers no special ad- 
vantage as a probe. We mention, however, the A c baryon, 
where the process A c — > t + K° allows for r in final state. 



The process in principle probes cus 



interaction. 



However, experimental study poses a problem. Not only 
would A c be hadronically produced, having a r in the 
final state does not allow for full reconstruction. For this 
reason, we shall drop r from consideration even if it is 
kinematically allowed in decay final state. For a rough 
estimate of the branching ratios in Eqs. i|36|) . (|37|l and 
we use for C^ cs£ and C^ cde respectively C^ UST and 
Cuudr f rom Eqs. dJl and (|rU|) obtaining 

B(D+ -► ~A£+) < 1.3 x 1CT 29 , (39) 
B(D° -> S+^+) < 1.0 x 10" 29 , (40) 
B{D°^p£ + ) < 4.0 xlO" 39 . (41) 

We now turn to the B system. Having dropped r from 
the final state, the transitions of interest are b — > cc£, 
cut, uu£, where I = /i, e. In a way, b — > \xcc is analogous 
to b — > sss, as far as generation counting is concerned. 
This provides some extra motivation, as CP violation in 
B — ► 4>Ks could be hinting at presence of new physics. 
One can in principle have the double charmed baryon 
S cc in final state, i.e. B°' + — > The S+ state is 

seen [2lJ in A+K~tt + channel, although it still needs to 
be confirmed. By analogy with D + — > K~tt + tt + which 
has ~ 9% rate, we expect #(S+ — > A+K~n + ) could be of 
order 20%. Since the cc system is generated from b — *• cc£ 
transition with nonrelativistic motion, there should be 
little suppression in the B — * S cc form factor. Again, 
we make analogy with r — > A7r~, D — > A£ + couplings of 
Eqs. © and l|35|) . and estimate the branching ratio as, 



The equivalent «b^h cc factor is set to a p value. Using 
Eq. we obtain 



r (5) 



ccbfj, 



, ca^ ; 8 !""^ < IP' 12 TeV 2 , (43) 

and from Eqs. (|4*2l and JlSJl we can estimate B(B°> + — > 
< 1.0 x 10- 28 . 

Although S cc decay is likely dominated by only a few 
channels, its reconstruction is still a problem. Typical A c 
decay modes are a few percent, so one would be limited 
by A c reconstruction efficiency. Perhaps one can also 
consider double charm (AC = 2) final states such as 
B + — * D* + D°p£~ , which is quite close to threshold, but 
would be stunning if seen: proton, £~ and (slow) 7r + from 
B vertex, while one has two D, rather than DD, tags. 
One could also consider the modes B — > A+ 1~ . 
Again, it would be hard to set a bound on the effective 
coupling, but can be searched for. 

For b — * cu£ generated processes such as B° — > A+^ _ , 
estimation becomes unreliable because of relatively large 
energy release [23, and the A c would be in relativistic 
motion in the B decay frame. But if we continue with 
our approach, the branching ratio is, 



B{B° -> A+T) < 1.1 x 10" 



\C UC M | 2 , 



(44) 



where a^^A c < ot p is applied. The suppression from 
a B^A c l a p cou ld be as large as 100, if one compares the 
observed B° -> A c p rate H3 with that of D+ir- . This 
would give B(B° — > A+£~) <4x 10 -30 , if one applies 
the bound of Eq. lfT5|). 

What may be more reasonable for the -B system, in 
view of the large available energy, is to adopt an inclusive 
approach. Note that b — ► cu£ _ bears some similarity 
with the usual b — > c£ + v, except for the wrong charge 
combination for ct~ vs c£ + arising from b decay (plus 
replacing v by u quark). With 6-tagging, one searches 
for B — > A+£rX, with no missing energy. If this can be 
properly achieved, it can be equated with the inclusive 
b — > cu£~ branching ratio, 



6(6 -> cut') = 4.0 x 10^ 2 |C.^ 



(45) 



and derive a more stringent bound. Applying the bound 
of Eq. (T5J}, we find B(b -> cu^~) < 1.3 x 10~ 27 . We re- 
mark that for b — > cu€ _ , because of unequal flavor index, 
the O 3 operator of Eq. © contributes. 

For inclusive 6 — > uu£~ transitions, we get, 



B(b^uu£-) = 7.4 x 10" 



\Cuube \ 2 - 



(46) 



B(B°> 



-r)<i.oxio- 4 |cf cM | 2 . 



(42) 



Exclusive search for _B — * p£~ can be done, but would 
not be profitable as this mode would be highly form fac- 
tor suppressed. This is because of the large energy re- 
lease [23 an d highly relativistic proton motion. However, 
what one should search for is B — > p£~X with 6-tagging, 
where p£~ has the wrong charge combination and come 
from the B decay vertex. Charm veto for semileptonic 



5 



decays should be applied. The study would probably be 
background limited due to mistag and vertex resolutions, 
but certainly should be pursued. The limit from Eq. , 
after applying the bound of Eq. (|Fo)l . is 2.4 x 10~ 27 . 

VI. CONCLUSION 

We have used experimental data on nucleon decays 
to give an estimate of BNV couplings involving higher 
generation leptons and quarks. Using an approach very 
similar to the one proposed by Marciano to infer very 
small branching ratio for BNV r decays, we find that 
BNV processes involving higher generation quarks are 
also extremely suppressed. 

Despite our findings, we believe it is still worthy to 



look for BNV processes in r, charm, B, and maybe in 
the future in top decays. Although we find that all BNV 
processes involving r and higher generation quarks are 
too strongly suppressed by proton stability, such that 
there seems no hope for observation at any future ma- 
chine, redundancy is very important. In this vein, we 
note that the CLEO Collaboration performed t~ — > jot , 
prj, p7r°7r°, pTT°T] : pj |13| search after the remark made 
by Marciano. The search has been extended recently to 
r~ — > A(A)7r~ by the Belle Collaboration. 
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